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The impact of nuclear physics theories on cooling of isolated neutron stars is analyzed. 
Physical properties of neutron star matter important for cooling are reviewed such as 
composition, the equation of state, superfluidity of various baryon species, neutrino emis- 
sion mechanisms. Theoretical results are compared with observations of thermal radiation 
from neutron stars. Current constraints on theoretical models of dense matter, derived 
from such a comparison, are formulated. 

1. INTRODUCTION 

Our knowledge of neutron star (NS) interiors is currently uncertain. In particular, 
the fundamental problem of the equation of state (EOS) at supranuclear densities in NS 
cores is still unsolved. Microscopic theories of dense matter are model dependent and 
give a large scatter of possible EOSs (e.g., Ref. jjj), from stiff to soft ones, with different 
compositions of inner NS cores (nucleons, hyperons, pion or kaon condensates, quarks). 
Thermal evolution of NSs depends on a model of dense matter which enables one to 
constrain the fundamental properties of dense matter by comparing simulations of NS 
cooling with observations. 

NSs are born hot in supernova explosions, with the internal temperature T ~ 10 11 K. 
In about one minute after the birth a star becomes transparent for neutrinos generated 
in its interiors. In the following neutrino-transparent stage the star cools via neutrino 
emission from the entire stellar body and via heat transport through the envelope to the 
surface and subsequent thermal surface emission of photons. In a hundred of years the 
NS crust and core become thermally adjust and the NS interior becomes isothermal (with 
the only temperature gradient located near the surface). After that the effective surface 
temperature, T s , reflects the thermal state of the core. 

The recent development of the theory has been reviewed, e.g., in Refs. Here, we 

outline the current status of the problem. 



2. OBSERVATIONS 

Observations of isolated NSs, whose thermal surface radiation has been detected or 
constrained, are summarized in Fig. ^ (following Ref. 0]). We present the estimated NS 
ages t and effective surface temperatures T s °° (as detected by a distant observer). 
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Figure 1. Observational limits of sur- 
face temperatures for several isolated NSs 
compared with the basic theoretical cooling 
curve of a non-superfluid NS model. 



Figure 2. Internal and surface tempera- 
tures; neutrino, photon and total luminosi- 
ties (redishifted for a distant observer) for 
the same NS model as in Fig. ^ 



For the two youngest sources only upper limits on the surface temperature T s °° have 
been established p, |6j. The surface temperatures of the next five sources, with ages 
10 3 ^ t < 10 5 years, have been obtained 0,0,0? U| by fitting their thermal radiation 



spectra with hydrogen atmosphere models. Such models are more consistent with other 
information on these sources (e.g., Ref. 0|) than the blackbod y m odel. On the contrary, 
for Geminga and PSR B1055-52 we present the values of T s °° (ill . 12T ] inferred using the 
blackbody spectrum because this spectrum is more consistent for these sources. The 
surface temperature of RX J1856. 4-3754 is still uncertain. Following |4| we adopt the 
upper limit T s °° < 0.65 MK. Finally, T s °° for RX J0720.4-3125 is taken from Ref. [3, 
where the observed spectrum is interpreted with a model of a hydrogen atmosphere of 
finite depth. 

As seen from Fig. ^ observational limits scatter in the — t plane. What can be 
learnt on dense matter in NS interiors from this scatter? 



3. THEORY VERSUS OBSERVATIONS 

A neutron star consists of a thin crust (of mass < 10 _2 M Q , where M Q is the solar 
mass) and a core (e.g., Ref. Q). The core-crust interface is placed at the mass density 
p ~ Po/2, where po ~ 2.8 x 10 14 g cm" 3 is the density of saturated nuclear matter. 
The crustal matter contains atomic nuclei, electrons, and (at p > 4 x 10 11 g cm" 3 ) free 
neutrons. The core is further subdivided into the outer (p < 2po) and inner parts. The 
outer core consists of neutrons, with an admixture of protons, electrons, and muons. The 
composition of the inner core is still unknown. It may be the same composition as the 
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Figure 3. Left: Illustrative models of critical temperatures for proton (p) and neutron (nt) 
pairing in NS core. Right: Neutrino emissivity in the same NS core at the temperature 
T = 3 x 10 s K for non-superfluid matter (thick line; noSF) and in the presence of either 
proton pairing (p) or proton and neutron pairing (p+nt). The vertical dotted line indicates 
the threshold of the direct Urea process. 



outer core but may also contain hyperons, pion or kaon condensates, quark matter, or a 
mixture of different phases. 

NS matter is strongly degenerate. The EOS, NS masses M and radii are almost temper- 
ature independent. Low-mass NSs (M ~ M Q ) have rather low central densities p c <2pQ 
and do not possess inner cores. NSs with masses close to the maximum allowable mass 
(M max ~ (1.5 — 2.5) Mq, for different model EOSs) have massive inner cores. 

NS cooling is calculated with a cooling code (e.g., Ref. in the form of cooling 

curves, T£°(t) (e.g., Fig. GJ. The initial cooling stage, t < 100 years, is accompanied by 
thermal relaxation of NS interiors (Fig. EJ). As long as t < 10 5 years, a star cools mainly 
via neutrino emission from its interiors (mainly from the core); this is the neutrino cooling 
stage. Later, at t > 10 5 years, the neutrino emission becomes inefficient, and the star 
cools via thermal surface emission of photons (the photon cooling stage). 

NSs may have different masses, surface magnetic fields, composition of surface layers, 
etc., but they are supposed to have the same EOS and superfluid properties of internal 
layers. In the absence of exact microscopic theory of NS matter we will use several model 
EOSs and phenomenological superfluidity models. 

The main regulators of NS cooling are: 

(a) EOS and composition of NS cores which affect neutrino emission mechanisms; 

(b) Superfluidity of baryons in NSs — it affects neutrino emission and heat content; 

(c) The presence of light elements (accreted envelopes) and strong magnetic fields in NS 
surface layers. These factors affect the thermal conductivity and the relation between the 
internal and surface temperatures of the star. 
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We discuss these regulators below (except for magnetic fields whose effects are examined 
in [Hj]). Other regulators are reviewed, e.g., in Ref. 0. 

Figure [T] shows the basic cooling curve. It is calculated for a star with a non-superfluid 
nucleon core, where the powerful direct Urea process of neutrino emission is forbidden. 
Such a star cools mainly via neutrino emission produced by the less powerful modified 
Urea process; the accretion envelope is absent. This basic curve is universal, being almost 
independent of the EOS and NS mass. It cannot explain all the observations - some NSs 
are hotter and some colder than predicted by the curve. However, one can explain the 
data by employing other cooling regulators. 

The effects of superfluidity and the direct Urea process, are demonstrated in Fig. El Here 
we adopt a moderately stiff EOS of dense nucleon matter suggested in Ref. [3| (the same 
version as used in 0). This EOS opens the direct Urea process at p > p D = 7.851 x 10 14 
g cm -3 , i.e., at M > = 1.358 M (M being the gravitational mass) and gives NS 
models with M max = 1.977 M Q . In non-superfluid matter the direct Urea process switches 
on sharply at p > pjy. However, neutrons and protons (like other baryons) in NS cores can 
be in superfluid state. As a rule, neutrons undergo triplet-state pairing, whereas protons 
undergo singlet-state pairing (e.g., Ref. 0]) with density dependent entieal temperatures 
^cnt(p) an d T cp (p) which are extremely sensitive to theoretical models. Superfluidity sup- 
presses traditional neutrino emission mechanisms (the modified and direct Urea processes 
and nucleon-nucleon bremsstrahlung) but opens a new neutrino process associated with 
Cooper pairing of baryons [l^- Figure El shows some phenomenological T cnt (p) and T cp (p) 
curves (from Ref. 0) and demonstrates that the direct Urea process and superfluidity 
greatly affect the neutrino emission (and, hence, NS cooling, as discussed later). 

The cooling can be strongly different for low mass, medium mass, and high mass NSs. 



3.1. Cooling of low-mass stars 

Low-mass NSs possess only outer nucleon cores. Some cooling curves are presented in 
Fig. 0] for NS models constructed with the same EOS as in Fig- El 

The solid curves are calculated assuming strong proton superfluidity p. It suppresses the 
modified Urea process in a low- mass NS. The neutrino luminosity of the star becomes lower 
(Fig. EJ), being determined by a weaker mechanism of neutrino emission (neutron-neutron 
bremsstrahlung, unaffected by superfluidity as long as neutrons are non-superfluid). This 
rises the cooling curves at the neutrino cooling stage. The thick solid curve is calculated 
for a star without any accreted envelope. This curve (contrary to the basic curve in Fig.[T} 
goes high enough to explain the sources hottest for their age (RX J0822-43, IE 1207-52, 
PSR B1055-52). Thus, we may treat these sources as low- mass NSs. The thin solid curve 
is calculated assuming, additionally, the presence of hydrogen or helium accreted envelope 
of mass AM = 1O -8 M . Light elements increase the thermal conductivity of the surface 
layers which further rises T£° at the neutrino cooling stage. The thin solid curve is close 
to the highest cooling curve provided by the standard cooling theory of NSs. 

Actually, we do not need very strong proton superfluidity (such as model p) to interpret 
the observations. The dashed lines in Fig. 0] are the same as the solid lines, but the proton 
critical temperature T cp (p) is reduced by a factor of 6. Weaker proton superfluidity relaxes 
the superfluid suppression of the modified Urea processes and lowers the cooling curves 
(with respect to the solid ones). Nevertheless, such superfluidity is still sufficient to 
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Figure 4. Cooling curves for the star with 
M = 1.3 M & and nucleon core compared 
with the observations (see text for other ex- 
planations) . 



Figure 5. Temperature versus depth z in 
the NS envelope for two internal temper- 
atures: 10 8 - 5 K (light curves) and 10 7 ' 5 K 
(heavy curves), for nonaccreted and ac- 
creted envelopes. 



interpret the observations: old hot sources are consistent with the thick dashed line (no 
accretion envelope) while young hot sources are well explained assuming an accreted 
envelope (thin dashed line). Finally, the dotted curve is the same as the thin dashed 
curve but the mass of the accreted envelope is assumed to decrease with time (e.g., due 
to diffuse burning 0) as AM(t) = AM(0) exp(-t/r), with r = 4000 years. This model 
can also explain the observations of NSs hottest for their ages. 

The effect of light-element envelopes on the cooling can be understood from Fig. It 
shows the growth of temperature within the NS of mass M = 1.4 M & and radius R = 10 
km. The solid lines refer to a non-accreted (Fe) surface. The dot-and-dashed lines are 
for the thickest accreted (He and C) envelope that can survive with respect to nuclear 
burning. Because He and C have higher thermal conductivity than Fe, their presence 
makes the NS envelope more heat-transparent, increasing T s . 

3.2. Cooling of high-mass neutron stars 

High-mass NSs have large central densities, masses M ~ M max , and contain massive 
inner cores. Microscopic theories predict that, as a rule, superfluidity dies out (and does 
not suppress neutrino emission) in the central parts of such stars. 

One can propose very different cooling scenarios of high-mass NSs. The simplest sce- 
nario assumes non-superfluid nucleon cores where the direct Urea process is forbidden. 
The corresponding cooling curves would be the same as the basic curve in Fig. ^ they 
cannot explain all the observations. 
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Figure 6. Left: A sketch of the neutrino luminosity L v versus stellar mass for NSs with 
the internal temperature T = 3 x 10 8 K at four models of NS structure. Right: Four 
hatched regions of T s °° which can be explained by cooling of NSs of different masses for 
four models in the left panel. 



It is widely thought that the neutrino emission in high-mass NSs is enhanced as com- 
pared to the emission provided by the modified Urea process. An enhanced emission 
would lead to fast cooling, allowing one to explain the observations of NSs coldest for 
their ages. There are different enhancement levels for different models of NS internal 
structure. Four scenarios are presented in Fig. El The left panel is a rough sketch of the 
neutrino luminosity as a function of M at T = 3 x 10 8 K. One can generally assume a 
slow neutrino emission in low- mass NSs (e.g., provided by the neutrino bremsstrahlung 
in neutron-neutron collisions), an enhanced neutrino emission in high-mass NSs, and the 
transition from the slow to enhanced emission with increa sing M in medium-mass NSs. 
The mass range of medium-mass stars is model dependent |2£|. In massive NSs L„ scales 
as T 6 for all scenarios except for Cooper pairing scenario (where L v oc T 8 0, 0|). In 
low- mass NSs, L v oc T 8 . 

The right panel of Fig.|U]sliows the limiting cooling curves for each scenario (no accreted 
envelopes). The upper curve refers to low- mass NSs. Their cooling is the same for all 
scenarios, and it explains the observations of NSs hottest for their age (Sect. 13. lj) . Four 
lower curves show cooling of maximum-mass NSs for four scenarios. They are the lowest 
cooling curves in these scenarios. The three lowest curves are taken from Ref. 0, while 
the fourth is from Ref. The range of T s °° between the upper curve and a lower curve 
can be filled by cooling curves of NSs with masses from ~ M Q to M max . Thus, we have 
four different ranges of T s °° (four hatched regions) for four scenarios. 

The highest enhancement of neutrino emission is provided by the direct Urea (Durca) 
process in nucleon (or nucleon-hyperon) cores |2j| . This scenario predicts the coldest mas- 
sive NSs and the widest theoretical T s °° range. If the direct Urea process is forbidden but 



7 



6.6 





lg t (yrs) 



Figure 7. Left: Cooling of NSs of several masses (indicated near the curves). NSs are 
assumed to have nucleon cores and proton superfluidity p from Fig. El Right: Same as in 
the left panel but adding the effect of neutron superfluidity nt. 



pion condensate is present in the inner NS core, the enhancement of neutrino emission is 
provided by the process of Durca type involving quasi-nucleons (e.g., 22\). This enhance- 
ment is weaker, the massive stars are hotter, and the acceptable range narrower. If 
pion condensate is absent, but kaon condensate available, the neutrino emission enhance- 
ment (in Durca- type processes involving quasi-baryons |22j) is even weaker and the T s °° 
range narrower. Nearly the same enhancement is expected in NSs with non-superfluid 
inner cores composed of quarks. Finally, the lowest enhancement can be produced in nu- 
cleon inner cores 0, 0| , where the direct Urea process is forbidden but mild superfluidity 
(e.g., of neutrons) is available (see Sect. 13.4)1 . It triggers the Cooper pairing neutrino 
emission which accelerates NS cooling. It gives the narrowest theoretical region of T s °°. 
As seen from Fig. |H1 all four scenarios are compatible with the observations. 



3.3. Cooling of medium-mass neutron stars 

The next question, crucial for explaining the observations (e.g., of the Vela and Geminga 
pulsars), is how cooling curves fill hatched regions in Fig. |U]if we vary the NS mass from 
~ M & to M max . The answer 23j is closely related to the contrast of slow and enhanced 
neutrino luminosities and the mass range of medium-mass NSs in Fig. El Let us outline 
this problem for NSs with nucleon cores. 

Figure [7| shows cooling of NSs of several masses with the same EOS as in Fig. |3] (no 
accretion envelopes). In the left panel we take into account strong proton pairing p, which 
extends to densities p > pr>. As long as T cp (p) > 3 x 10 9 K, it suppresses the modified 
and even the direct Urea process and leads to very slow neutrino emission. At higher p it 
gradually dies out, opening the direct Urea process. The gradual opening broadens the 
direct Urea threshold (Fig. |3J) and ensures the gradual decrease of cooling curves with 
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Figure 8. Left: Model density dependence of critical temperatures of protons (pi) and 
neutrons (ntl) in a nucleon NS core for the EOS which forbids the direct Urea process. 
Right: Cooling curves of NSs of several masses for the same EOS, taking into account 
superfluidities pi and ntl. After Ref. 0. 



increasing M. In this way we may attribute masses to observed NSs |24|. For instance, 
we obtain M « 1.47 M for the Vela pulsar. However, this weighing of NSs is sensitive to 
the EOS of dense matter, the threshold of the direct Urea process, and the superfluidity 
model T cp (p). Were superfluidity absent, the transition from slow to fast cooling would 
occur in a very narrow mass interval (0 < M — Md < 0.001 M ), and the successful 
interpretation of the data would be unlikely (e.g., Ref. [2J). 

3.4. Harmful and useful Cooper-pairing neutrino emission 

The cooling effect of neutrino emission due to Cooper pairing of nucleons may be 
different. For instance, let us assume the presence of neutron pairing nt with the peak 
of T cnt (p) as low as ~ 4 x 10 8 K at p ~ 4 x 10 14 g cm" 3 (Fig. EJ). This superfluidity is 
mild and insignificant, according to nuclear physics standards, but crucial for NS cooling. 
It appears in a cooling star when the internal temperature falls below the peak value. It 
creates then a powerful neutrino emission owing to Cooper pairing of neutrons in outer 
NS cores (especially efficient in low-mass NSs). The emission accelerates NS cooling, as 
shown in the right panel of Fig. and violates the interpretation of the observations 
of such sources as PSR B1055-52. Thus, this mild neutron superfluidity contradicts the 
observations. 

The opposite example is given in Fig. |HJ Let us consider NSs with nucleon cores and 
employ the EOS jl^] which forbids the direct Urea process in all NSs with M < M max = 
2.05 Mq. Furthermore, let us adopt the model of strong proton pairing pi and mild 
neutron pairing ntl (the left panel of Fig. |HJ). Pairing pi is similar to pairing p in Fig. 
El it suppresses the modified Urea process in low-mass NSs. The peak of T cnt (p) for 
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pairing ntl is as low as for pairing nt but shifted to higher densities. Accordingly, pairing 
ntl is inefficient in low-mass stars and does not speed up their cooling. However, the 
enhanced neutrino emission owing to this pairing operates in massive NSs and accelerates 
their cooling (a scenario considered in Sect. 13.21 Fig. Then, as seen in the right 
panel of Fig. |H1 the cooling of NSs of different masses enables us to explain the data. 
In this case mild neutron pairing is useful for interpretation of the observations but a 
successful interpretation is possible only under stringent constraints on the T cnt (p) profile 
Moreover, a discovery of a new NS slightly colder than those observed now would 
ruin this interpretation. 



4. CONCLUSIONS 

We have outlined several possible scenarios of NS cooling. In particular, we have con- 
sidered cooling of NSs with nucleon (nucleon/hyperon) cores, and cores containing exotic 
phases of dense matter. We have shown that many scenarios are currently compatible 
with observations of thermal radiation from isolated NSs. Our main conclusions are: 

(i) Some NSs (e.g., RX J0822-43 or PSR B1055-52) are hotter and some (e.g., the Vela 
pulsar) colder than non-superfluid NSs which cool via the modified Urea process. Hotter 
NSs are possibly low-mass stars, while coldest observed NSs are possibly more massive. 

(ii) Currently, the observations unable one to discriminate between many cooling sce- 
narios. However, they seem to rule out mild superfluidity with the peak of the critical 
temperature T c (p) between ~ 3x 10 8 and ~ 2 x 10 9 K at p < 8 x 10 14 g cm -3 in NS 
cores. This superfluidity would initiate Cooper pairing neutrino emission in low-mass 
NSs hampering the interpretation of the observations of old and warm NSs, such as PSR 
B1055-52. In contrast, mild superfluidity with the peak of T c (p) at higher p can be useful 
for interpretation of the observations. 

We have discussed main cooling scenarios but not all of them. Some others are reviewed 
in Refs. Actually, the effects of superfluidity are more sophisticated than discussed 

above. For instance, cooling curves do not change qualitatively by exchanging T c (p) for 
neutrons and protons [26j. Strong superfluidity of all baryon species (with peaks of T c (p) 
higher than 2 x 10 9 K) leads to a very low heat capacity of NSs. Such stars appear at 
the photon cooling stage earlier than at t ~ 10 5 years; they are too cold at that stage. 
Weak superfluidity, with T c (p) <3x 10 8 K, does not occur in NSs of ages t < 10 5 years 
and does not affect their cooling. Cooling of low-mass NSs can be strongly affected by 
singlet-state pairing of neutrons in inner stellar crusts. 

New observations of NSs are required for a better understanding of their internal struc- 
ture. New discoveries of cold NSs would be especially useful. Observations of cooling NSs 
can be analyzed together with other observational data, for instance, with observations 
of quiescent thermal emission from NSs in soft X-ray transients (see Q, for references). 
This would allow one to obtain more stringent constraints on NS structure. 
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